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We study the zero-temperature phase diagram and the low-lying excitations of a square-lattice spin- 
half Heisenberg antiferromagnet with two types of regularly distributed nearest-neighbour exchange 
bonds ( J > (antiferromagnetic) and — 00 < J' < 00) using the coupled cluster method (CCM) for 
high orders of approximation (up to LSUB8). We use a Neel model state as well as a helical model 
state as a starting point for the CCM calculations. We find a second-order transition from a phase 
with Neel order to a finite-gap quantum disordered phase for sufficiently large antiferromagnetic 
exchange constants J' > 0. For frustrating ferromagnetic couplings J' < we find indications that 
quantum fluctuations favour a first-order phase transition from the Neel order to a quantum helical 
state, by contrast with the corresponding second-order transition in the corresponding classical 
model. The results are compared to those of exact diagonalizations of finite systems (up to 32 sites) 
and those of spin-wave and variational calculations. The CCM results agree well with the exact 
diagonalization data over the whole range of the parameters. The special case of J' = 0, which is 
equivalent to the honeycomb lattice, is treated more closely. 



I. INTRODUCTION 

The subject of quantum spin-half antiferromagnetism 
in low-dimensional systems has attracted a great deal of 
interest in recent times in connection with the magnetic 
properties of the cuprate high-temperature superconduc- 
tors. However, low-dimensional quantum spin systems 
are of interest in their own right as examples of strongly 
interacting quantum many-body systems-. Although we 
know from the Mermin- Wagner theoremEJ that thermal 
fluctuations are strong enough to destroy magnetic long- 
range order at any finite temperature, the role of quan- 
tum fluctuations is less understood. As a result of inten- 
sive work in the late eighties, it is now well-established 
that the ground-state of the Heisenberg antiferromag- 
net on the square lattice with nearest-neighbour inter- 
actions is long-range ordered (see for example the re- 
view in Ref. However, Anderson's and Fazekas' 
investigationsa of the triangular lattice led to a conjec- 
ture that quantum fluctuations plus frustration may be 
sufficient to destroy the Neel-like long-range order in two 
dimensions. Another specific area of recent research is 
the spin-half J1-J2 antiferromagnet on the square-lattice 
where the frustrating diagonal Ji bonds plus quantum 
fluctuations are able to realize a second-order transition 
from Neel ordering to a disordered quantum spin liq- 
uid (see for example Refs. |4p|7] and references therein). 
On the other hand, there are cases in which frustration 
causes a first-order transition in quantum spin systems in 
contrast to a second-order transition in the correspond- 



ing classical model (see for example Refs. || p"T| ). 

Besides frustration, there is another mechanism to re- 
alize the "melting" of Neel ordering in the ground states 
of unfrustrated Heisenberg antifcrromagnets, namely the 
formation of local singlet pairs of two coupled spins. This 
mechanism may be relepirtLfor the quantum disordered 
state in bilayer systemsE3~li3 as well as in CaV^g (see 
for example Refs. |5~6| — p~S| and references therein). The 
formation of local singlets is connected with a gap in the 
excitation spectrum. By contrast, the opening of a gap 
in the excitation spectrum of frustrated systems seems 
to be less clear and might be dependent on details of the 
exchange interactions. 

In the present paper, we study a model which contains 
both mechanisms, frustration and singlet formation, in 
different parameter regions. We mainly use in this arti- 
cle the coupled cluster methodJijiEj which has become 
widely recognized as one of the most powerful and most 
universal techniques in quantum many-body theory. In 
recent years there has been increasing apcess in apply- 
ing the CCM to quantum spin systems, Q'E3~E3 especially 
with the advent of high-prder approximations which uti- 
lize computer algebra.Ba Subsequently, high-order CCM 
approximations have been applied to the XX Z modgl,E2l 
the anisotropic XY modeled and the J1-J2 model.0 In 
addition to the CCM results we also present variational, 
spin-wave theory (SWT) and exact diagonalization (ED) 
results for the sake of comparison. 
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II. THE MODEL 



We consider a spin-half Hcisenbcrg model on a square 
lattice with nearest neighbor bonds J and J' in a regular 
zigzag pattern as shown in Fig. [ll The Hamiltonian is 
given by 



J=i 



H = J Si ■ Sj + J' 2_. S i ' S j 

<ij>l <ij>2 

EE( 1 + ^( J '- 1 )) s '' s ^- w 

The sums over < ij > 1; and < ij > 2 represent sums over 
the nearest-neighbour bonds shown in Fig. |l| by dashed 
and solid lines, respectively. Throughout the paper we fix 
the J bond to be antiferromagnetic (J > 0) and hence- 
forth scale it to the value J = 1, and consider J' as the 
free parameter of the model. We also split the square 
lattice into the equivalent A and B sublattices shown in 
Fig. |]. In Eq. (Q) the sum over i runs over the sites of 
the sublattice A, with vectors p — {(0, ±1), (±1, 0)} con- 
necting nearest neighbours. In particular, pj> — (1,0) 
pertains to the coupling with J' bonds. 

Each square lattice plaquette consists of three J = 1 
bonds and one J' bond. In the case of ferromagnetic 
J' bonds (i.e., J' < 0), the plaquettes are frustrated. 
Conversely, for antiferromagnetic J' bonds (i.e., J' > 0) 
there is no frustration in the system, although the differ- 
ence of the coupling strengths J and J' leads to quantum 
competition. This modeLlias been treated previously us- 
ing perturbation theory,Ld renormalized spip_wave theory 
(RSWT)El and exact diagonalization (ED) £3 It allows us 
to study the influence of local singlet formation (J' > 1) 
and frustration (J' < 0) on the stability of the Neel order 
within a single model. 

Ferromagnetic bonds in an antifectoaiagnetic matrix 
have been discussed in recent timesQ'cJtJ iiir-epnnection 
with the proposal by Aharony and coworkersell to model 
localized oxygen holes in the Cu-O-planes by local ferro- 
magnetic bonds between the copper spins. It was argued 
that random ferromagnetic bonds may influence the an- 
tiferromagnetic order drastically and ma$ -support the 
realization of a quantum spin-liquid stateO'EErEfl 

On the other hand, the case of antiferromagnetic J' 
bonds with J' > 1 resembles the situation in bilayer sys- 
tems and in the depleted square-lattice antiferromagnet 
CaV40g, in which the competition between two different 
antiferromagnetic bonds leads to a phase transition from 
antiferromagnetic long-range order to quantum disorder 
with a finite gap. It is seen in this article that the transi- 
tion point obtained for the model of Eq-ipl) is quite close 
to that obtained for the bilayer model.Eij 

There are some special cases of the model Hamiltonian 
ofEq. (|): 

(i) J' = 1: square-lattice antiferromagnet, for which 
the ground state is long-range ordered; 

(ii) J' = 0: honeycomb-lattice antiferromagnet, for 



which the ground state is long-range ordered; 

(iii) J' = — oo: spin-1 triangular lattice, for which the 
ground state is long-range ordered; and 

(iv) J' = +oo: valence-bond solid, for which the ground 
state is a rotationally invariant quantum dimer state with 
an excitation gap. 

Classical ground state. For J' > —1/3 the Neel state is 
the classical ground state of the Hamiltonian of Eq. (Q). 
At J' c — —1/3 there is classically a second-order phase 
transition to a ground state of helical nature (see Fig. [I]), 
with a characteristic pitch angle $ = ±|$ c i| given by 



(iA 



J' > 
j' < 



(2) 



where tip, different signs correspond to the two 
chiraliticsE3 of this helical state. Note that for $ = 
this is just the Neel state. More generally, the pitch 
angle varies with J' from |$ c i| = for J' > —1/3 to 
|$ cl | = vr/3 for J' = -oo. Note that |$ cl | = vr/3 (re- 
alized at J' = —oo) corresponds to the ground state of 
the spin-1 triangular lattice. We describe the directions 
of the spins and belonging to the A and B sub- 
lattices respectively, for the classical helical state with a 
characteristic angle $ as followsEll (and see Fig. ([!])), 

(R) = u cos Q • R + v sin Q • R, (3) 
s B (R + x) = ucos(Q • R + vr + 3$) + vsin(Q ■ R + n + 3$), 

where u and v are perpendicular unit vectors in the spin 
space, R runs over the sites of the sublattice A, and we 
have Q = (2$, 0) for the pitch vector Q. We note that 
this general helical state does not have a periodicity in the 
x-direction because $ is in general not of the form rmr/n 
with m and n integral. We also note that we have only 
three different angles between nearest-neighbour spins, 
namely ±(tt — $) for the J = 1 couplings and it — 3$ for 
the coupling with J'. 

The maximum frustration is in the region around 
J' ps —1. Bearing in mind the situation for the J\-Ji 
model, one might expect that for the extreme quantum 
case (spin-half) quantum fluctuations might be able to 
open the window to a spin-liquid phase for a finite range 
of parameters around this region of maximum frustra- 
tion. On the other hand, for strong antiferromagnetic 
bonds (J' ^> 1) there is, of course, no indication in the 
classical model for the breakdown of the Neel order. 

Simple variational ansatz for the quantum ground- 
state. In the quantum case, the region of strong anti- 
ferromagnetic J' bonds (J' 3> 1) is characterized by a 
tendency to singlet pairing of the two spins correspond-, 
ing to a J' bond. Using a high-order series expansionc2l 
the Neel order was found to be stable up to a critical 
value J' s w 2.56. 

A comparable value can be obtained using j-a simple 
variational wave function similar to that usedtil for bi- 
layer systems, namely 
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[I ^4.*+*)-*] 4.iT*+4>], (4) 



where the lattice sites i and i + x correspond to a J' 
bond, and where the product in Eq. (|J) is thus effectively 
taken over the J' bonds of the lattice of Eq. ([[]). The 
trial function depends on the variational parameter t and 
interpolates between a valence-bond state realized for t = 
1 and the Neel state for t — 0. For t = 1, the singlet 
pairing is complete and l^var) represents an eigenstate 
of the model of Eq. (Q) in the limit J' — * oo (dimer 
state). By minimizing (^varl^l^var) with respect to the 
variational parameter t we get an upper bound for the 
ground-state energy per spin of the model of Eq. (|l|) , 



-(J' 2 + 3J' + 9)/24 
-3J78 



J' < 3 
J' > 3 



(5) 



The relevant order parameter describing the Neel order 
is 



M =(* 1^1* ) _{l/Vl-^' 2 /9 J'<3 

Ju var — \*var|*i l^var/ — \ q J' > 3 ' 

(6) 

showing a breakdown of the Neel order at a critical value 
J' = 3. 



III. COUPLED CLUSTER CALCULATIONS 
A. The ground-state formalism 

The starting point for any CCM calculation (see over- 
wiev in Ref. |^) is the choice of a normalized model or 
reference state together with a set of mutually com- 
muting multispin creation operators CJ which are de- 
fined over a complete set of many-body configurations 
/. The operators Cj are the multispin destruction op- 
erators and are defined to be the Hermitian adjoints of 
the CJ. We choose {|$);C/} in such a way that we 
have ($|C/ = = Cj|$), V/ ^ 0, where, by definiton, 
Cq = 1, the identity operator. 

For spin systems, an appropriate choice for_ the CCM 
model state |$) is often a classical spin statejEl in which 
the most general situation is that each spin can point in 
an arbitrary direction. For the case of the Hamiltonian of 
Eq. (0), we choose the helical state illustrated in Fig. [TJto 
be our model state. Although the classical ground state 
of Eq. (0) is precisely of this form, we do not choose the 
classical result for the pitch angle <& but we consider it 
rather as a free parameter in the CCM calculation. 

In order to perform a CCM calculation, we would like 
to treat each site equivalently and we do this by perform- 
ing a rotation of the local spin axes at each site about 
the y-axis such that all spins in the model state align in 
the same direction, say down (along the negative z-axis) . 
After this transformation wc have 



(<) 



(where the indices r, r', r" , . . . denote any lattice site) re- 
spectively, for the model state and the multispin creation 
operators, which now consist of spin-raising operators 
only. 

In order to make the spin Sj point down let us suppose 
we need to perform such a rotation of the spin axes by 
an angle Si. This is equivalent to the transformation 



cos 

of 



SiS* + suites- 



sf — > — sin Sisf + cosSisf. 



(8) 



A similar rotation about the y-axis by an angle Sj is 
performed for the spin Sj. Thus we get for the transfor- 
mation of the scalar product of the two spins, Sj • Sj — > 
(si -Sj) v , where 

(si • Hj) v = sin - af sf] + cos^[s?s| + sis'] + aftf 



+j(cos(p + l)[sfs j +s 



4 

+- (cos <p - l)[sf 
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sj + 8i 8j 



COS (yOSj Sj 

(9) 



The angle (p = Sj — Si is the angle between the two spins, 
and s ' = s x ± is v are the spin-raising and spin-lowering 
operators. Note that this product of two spins after the 
rotation depends not only on the angle between them, but 
also on the sign of this angle. In case of the Neel model 
state ($ = 0), the angle between any neighbouring spins 
is 7r, and hence Eq. (m becomes Sj • Sj - 



2 l b i b j 



Using the helical state of Eq. (|3j) with the characteristic 
angle the Hamiltonian of Eq. (|l]) is now rewritten in 
the local coordinates as, 



(10) 



i£A p 



where the angles between neighbouring spins are ip±y = 
7r + $, ip-£ = 7T — $ and (fix — 7r + 3$. While the general 
helical state (see Fig. [I]) does not have translational sym- 
metryin the x-direction, the transformed Hamiltonian of 
Eq. (|1C|) does have this symmetry since it depends only 
on the angles between neighbouring spins. 

Having defined an appropriate model state |$) with 
creation operators Cf , the CCM parameterizations of 
the ket and bra ground states are given by 



I*) =e S |$), 



S = ^S I C+, 

1*0 



S = l 



^s lCl . 

1*0 



(11) 



(12) 
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The correlation operator S is expressed in terms of the 
creation operators CJ and the ket-state correlation co- 
efficients Si. We can now write the ground-state energy 
as, 

E = (<l>\e- s He s \<i>). (13) 

To describe the magnetic order of the system, we use a 
simple order parameter which is expressed in terms of 
the local, rotated spin axes, and which is given by 

M = -<*|*f|*>, (14) 

such that the order parameter represents the on-site mag- 
netization. Note that M is the usual sublattice magneti- 
zation for the case of the Neel state as the CCM model 
state. 

To find the ket-state and bra-state correlation coef- 
ficients we have to require that the expectation value 
H = (i&\H |\&) is a minimum with respect to the bra-state 
and ket-state correlation coefficients. This formalism is 
exact if we include all possible multispin configurations 
in the correlation operators S and S, which is usually 
impossible in a practical situation. We use the LSUB?i 
approximation schemecS to truncate the expansion of S 
and S in the Eqs. (|ll|) and 

Using the lattice symmetries, we have now to find 
all different possible configurations with respect to the 
point and space group symmetries of both the lattice 
and Hamiltonian with up to n spins spanning a range of 
no more than n adjacent lattice sites (LSUBn approx- 
imation) and these are referred to as the fundamental 
configurations. 

The Hamiltonian of Eq. (|l|) has four lattice point-group 
symmetries namely two rotational operations (0°, 180°) 
and two reflections (along the x- and y-axes), defined by: 

x^x, y^y; x -> -(x + 1), y -> -y, 
x^x, y^-y; x^-(x + l), y -> y. (15) 

The rotation of 180° and the reflection along the y-axis 
are connected by a shift of x — (1,0). The translational 
operator T is defined by 

T = (n + m)x + (to — n)y, n, m integral, (16) 

such that translational symmetry is preserved. 

The Neel model state also contains these symmetries, 
and so for this model state we can directly apply all these 
symmetries in finding the fundamental configurations. 
On the other hand the general helical model state ($ ^ 0) 
has only two of the above four lattice point-group sym- 
metries, namely x — > x, y — > y and x — ► x, y — ► —y, 
and so this reduced symmetry yields a larger number of 
fundamental configurations. 

In the case of the Neel model state (<J> = 0), the num- 
ber of fundamental configurations can further be reduced 
by explicitly conserving the total uniform magnetization 



s?p = J2k s t (* ne sum on ^ runs over au lattice sites) be- 
cause the ground state is known to lie in the = sub- 
space. This means that we exclude configurations with an 
odd number of spins, and therefore we do not use LSUB3, 
LSUB5, etc. approximations. The helical state is not an 
eigenstate of Sj, and we cannot apply this property when 
using the helical model state. The fundamentaLconfigu- 
rations can now be calculated computationally,!^! and the 
resulting numbers of LSUBn configurations for n < 8 are 
given in Table |. 

The Ket-stater-ajid bra-state equations are calculated 
computationally.!^! For the Neel model state, we are able 
to carry out the CCM up to the LSUB8 level (where we 
need to solve 4986 coupled equations), whereas for the 
helical state we could do this only up to the LSUB6 level 
(where we need to solve 1638 coupled equations). 

B. The excited state formalism 

We use the excited-state formalism of Emrich@IE3'E3 
to approximate the excited-state wavef unctions. We ap- 
ply an excitation operator X e linearly to the ket state 
wavefunction (|TT|), such that 

|* e ) = A e e s |$); X e =J2 X I°t- ( 17 ) 

Using the Schrodinger equation, H\^> e ) = I? e |^ e ), we 
find that 

e e X e |$) = e- s [H, X e ]_e 5 |$), (18) 

where e e (= E e —E) is the difference between the excited- 
state energy (E e ) and the ground-state energy (E). Ap- 
plying ($|C/ to Eq. dH) we find that, 

e e X? = ($|C /e - s [ff,X e ]_e 5 |$), (19) 

which is an eigenvalue equation with eigenvalues e e and 
corresponding eigenvectors Xf . 

As for the ground state, we must use an approximation, 
scheme for X e in Eq. ( p7| ) . Although it is not necessaryc3 
to use the same approximation for the excited state as 
for the ground state, we in fact do so to keep the CCM 
calculations as systematic and self-consistent as possible. 
We define the fundamental configurations for LSUBn (for 
the Neel state) as previously, but we now restrict the 
choice of configurations to contain only those which pro- 
duce £u change of sf, of ±1 with respect to the model 
state.ufj Since we are only interested in the lowest-lying 
excitations, the restriction to these single-magnon spin- 
wave-like excitations is the correct choice. The number 
of fundamental excited-state configurations for LSUBrt is 
given in Table ffl. 

To calculate the terms of the right hand side of Eq. ( |l9| ) 
we use the same computational algorithm as for the cal- 
culation of the ground-state ket equation. The terms 
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contain the ground ket-state correlation coefficients Si, 
so once these coefficients have been determined the eigen- 
value equation ([l9]) can be solved (numerically) . We fur- 
thermore choose the lowest energy eigenvalue of Eq. (|l9| ) 
in order to calculate the excitation energy gap, A. We 
note that the eigenvalues of Eq. (|l9|) are not guaranteed 
to be real, since as a generalized eigenvalue equation it is 
not symmetric. However, over the entire regime of inter- 
est, the values of A so obtained are found to be real. We 
have performed these calculations for the excited state 
up to the LSUB6 level of approximation. 



C. Extrapolation of the CCM-LSUBti results 



Although no scaling theory for results of LSUBn ap- 
proximations, ha s.. yet been proven, there are empirical 
indicationsEJn3'E£rc3 of scaling laws for the energy, the 
magnetization, and the excited-state energy gap for var- 
ious spin models. 

These scaling laws can be justified by the observa- 
tions that they fit the results well (i.e., with low mean- 
square deviation), and that the extrapolated results are 
in good agreement with results of other methods (e.g., 
Green function Monte Carlo or series expansion for the 
2D XJLZ modefiil'ES) or with exact results (e.g., ID XY 
modeled) . In accordance with those previous results we 
use the following scaling laws: for the ground-state en- 
ergy, 



(20) 
(21) 
(22) 



E = a + ai(l/n 2 ) + a 2 (l/n 2 ) 2 ; 
for the ground-state magnetization, 

M = b + b 1 (l/n) + b 2 (l/n) 2 ; 
and for the gap of the lowest-lying excitations, 

A = c + Cl (l/n) + c 2 (l/n) 2 ; 
where n is the LSUBn approximation level. 

D. Choice of the CCM model state 



As stated previously, we use the helical state of Eq. (||) 
with the characteristic angle <£>, illustrated in Fig. [j, as 
the model state for the CCM. We must therefore make a 
selection of an appropriate value for $. A possible choice 
would be the classical ground state of the Hamiltonian 
of (0) (i.e., $ = $ c i as given by Eq. @)). 

Another possibility is to perform a CCM-LSUBn ap- 
proximation calculation and then to minimize the corre- 
sponding LSUBn approximation to the energy with re- 
spect to 



LSUB 



»(*) 



<^ $ = $ 



LSUBn ■ 



(23) 

The results for <&LSUBn will be given later (Fig. ^). How- 
ever, we note now that although the CCM does not yield 



a strict upper bound for the ground-state energy, using 
$ = ^LSUBn (he., using the CCM with a variational p*u 
rameter) has been found to be a reasonable assumption.^ 
There are several additional arguments to suggest that 
$ = ^LSUBn is indeed a better choice than $ = <I> c i, as 
indicated below. 

In the first place we note that we cannot find solutions 
for the LSUB6 equations using $ = <f> c i( J') in the region 
—0.7 < J' < —0.47 insofar as the Newton method used 
to solve these equations does not converge in that region. 
This is a clear indication that this model state is not a 
good one. By contrast, such behaviour is not found for 

* = *LSUBn- 

Secondly, it is generally known that quantum fluctua- 
tions tend to prefer collinear ordertre 2 ] (e.g., Neel order). 
We will indeed find (and see Fig. fj] below) that the Neel 
ordering ($ = 0) seems to survive for some J' < —1/3, 
in which region it has already broken down in the classi- 
cal case. This is also in agreement with results of exact 
diagonalizations for our model (and see Fig. || below). 

Thirdly we find better agreement of the CCM results 
for the energy compared to exact diagonalization results 
by using the helical state as the model state with the 
value $ = ^LSUBn rather than with the classical value 
<f> = <f>ci. We find that CCM results for the ground-state 
energy usually agree well with the corresponding ED re- 
sults (and with results of other methods), E3 provided that 
a good CCM model state is chosen. 

We therefore use the helical state with $ = ^lsub™ as 
the CCM model state throughout this paper. Note that 
for J' > —1/3 this model state is identical to the classical 
ground state of Eq. (fil) but that for J' < —1/3 it is not. 



IV. RESULTS 

Using the CCM scheme described above, we calcu- 
late the approximate ground state and the low-lying ex- 
citations of the Hamiltonian of Eq. (|l|). For compar- 
ison we also exactly diagonalize finite sized lattices of 
square shape. We use periodic boundary conditions with 
N = 16, 18, 20, 26 and 32 spins, and we extrapolate t^ha 
infinite system using standard finite-size scaling lawsO'O 
We present results for the ground-state energy, the order 
parameter and the excitation gap. We examine the for- 
mation of local singlets (for J' > 1), the effects of frustra- 
tion (for J' < 0), and the special case of the honeycomb 
lattice (J 7 = 0). 

J' > 1: Formation of local singlets. Using the CCM 
we obtain clear indications of a second-order phase tran- 
sition to a disordered dimer-like phase at a certain crit- 
ical value of J', namely J' s . For J' > J' s , the Neel- like 
long range order melts, (i^-, the sublattice magnetiza- 
tion, M, given by Eq. (|14j ) becomes zero). Our esti- 
mate for J' s using the four extrapolated LSUBn results 
for M with n = 2,4,6,8 (see Fig. |) is J' s « 3.41. How- 
ever, using only the three CCM LSUB71 approximations 
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with n = 4, 6, 8 for the extrapolation, we obtain a value 
J' s w 3.16, which indicates that the true value could be 
even somewhat smaller. This is in agreement with our 
corresponding result using exact diagonalizations of small 
systems. By using the extrapolation scheme of Ref. 
we find a critical value J' s 2.45 for the magnetization. 
Note however that better accuracy requires larger sys- 
tems because of the exact diagonalization (ED) extrapo- 
lation ansatz for M (i.e., M = + const x N^ 1 / 2 ). 
Therefore, we cannot consider the ED results for the 
magnetization (and see Fig. ^) as quantitatively correct. 
Our two results for the critical J' s also agree with the 
estimate J' s « 2.56 from series expansion,El and even 
the result J' s = 3 from the simple variational ansatz of 
Eq. (Q) agrees surprisingly well with these values. By 
contrast, the_aecond-order renormalized spin wave the- 
ory (RSWT)ELl gives the larger result J' s w 5.0, indicat- 
ing that the standard spin-wave approach is insufficient 
to describe this type of transition. 

Another indication of a dimcrized phase is the appear- 
ance of a gap A between the ground state and the lowest- 
lying excited state. We clearly expect a spectrum with 
gapless Goldstone modes if the ground state is Neel long- 
range ordered, whereas for a disordered singlet ground- 
state the formation of triplet excitations may cost a fi- 
nite amount of energy. This behaviour is reflected by our 
results using both CCM and exact diagonalization (see 
Fig. ||), which agree well with each other. For J' ^> J' s , 
there is a gap proportional to J', corresponding to the 
dimer-like nature of the ground state. The gap obvi- 
ously opens in the range 2.5 < J' s < 3.0 in both the ED 
and CCM calculations. This is in good agreement with 
the corresponding estimates for the critical point using 
the order parameter. Note, that the standard linear spin 
wave theory fails in calculating the gap ^e., gives gapless 
modes for all values of J' with J' > 0)M 

By comparing the results for the ground-state energy, 
we find excellent agreement between the CCM results 
and the results from exact diagonalization (see Fig. |J) 
for J' > Q_. By contrast, spin-wave theory (SWT) 
calculations^ show a significant deviation from these re- 
sults for larger J' >■ 1. These spin- wave results are obvi- 
ously poor since the simple upper bound for the energy 
given by Eq. (0) (e.g., Eq = —1.5 for J' — 4) is smaller 
than the corresponding SWT results (e.g., Eq = —1.42 
from 2nd-order RSWT). By contrast, CCM and ED re- 
sults (both are about Eq = —1.54 for J' = 4) are slightly 
smaller than the variational result. While both CCM 
and SWT calculations have the Neel state as starting 
point, we find the CCM is much better able than SWT 
to describe the transition to the rotationally invariant 
disordered state and to the completely dimerized state 
(represented by the variational function of Eq. (Q) with 
t = 1). Note that even the simplest CCM approxima- 
tion (LSUB2) gives the correct asymptotic result for the 
energy (i.e., Eq. (||)) for very large values of J', whereas 
SWT does not. 

For the case of the pure square-lattice Heisenberg anti- 



ferroma gnet (i .e., J' = 1), we reproduce the CCM results 
of Refs. £8|,fl0|, which have already been demonstrated to 
agree well with those from other methods. 

J' = 0: honeycomb lattice. For the special case of 
J' = (which is equivalent to the honeycomb lattice), 
we find that the CCM and the ED results are in good 
agreement (see Table However, the magnetization M 
for ED is found to be smaller then the CCM result. We 
note, however, that the CCM result for Jkf at this point 
agrees with the result of high-order SWTEJXM = 0.28) as 
well as with the result of series expansions (M = 0.26), 
although it does naLagree so well with the result of Monte 
Carlo calculationsQ (M = 0.22). We note too that our 
CCM results here agreej-oerfectly with previous lower- 
order CCM calculations.EZI 

J' < 0: Frustration. For J' < —2, we find that the 
extrapolated ED results for the energy lie appreciably 
above the CCM (and SWT) results (see Fig. @). This 
is because the energies for the small lattices considered 
do not fit well to the finite-size scaling law (E/N = 
Eoo/N + const x N~ 3 / 2 ) in this region. The finite-size 
effects for systems with an incommensurate helical struc- 
ture are found to be larger than for systems with, for 
example, Neel order or with dimerized spin pairs. How- 
ever, we find that our best ED result (with 32 spins) 
shows only very small deviations from the CCM result, 
even in the frustrated region. 

While classically we have a second-order phase tran- 
sition (from Neel order to helical order) at J' c = —1/3, 
using the CCM we find indications for a shift of this 
critical point to a value J' c ~ —1.35 (see Fig. [?]) in the 
quantum case. The ED data of the structure factors (see 
Fig. |B|) also show a shift of the transition to stronger fer- 
romagnetic J' bonds. Both correspond to the general 
picture that quantum fluctuations prefer a collinear or- 
dering (such as Neel order). Hence this ordered state 
can survive for the quantum case into a region where 
classically it is already unstable. The Neel model state 
($ = 0) gives the minimum ground-state energy for all 
values J' > J' c , where J' c is also dependent on the level 
of LSUBn approximation level. For J' < J' c another 
minimum in the energy for $ ^ is found to lie lower 
than the minimum at $ = (see Fig. ||). The state 
for $ ^ is believed to be a quantum analogue of the 
classical ground state in two dimensions. Furthermore, 
the crossover from one minimum solution to the other is 
not smooth but is abrupt at this point (see Fig. || and 
Fig. 0). This behaviour is assumed here to be an indi- 
cation of a phase transition. Furthermore, it might also 
indicate that this is a first-order phase transition and, 
consequently, that due to quantum fluctuations the na- 
ture of this phase transition is changed from the classical 
second-order type to a first-order type. 

The behaviour of the order parameter (see Fig. ||) 
in the region around J' c ~ —1.35 where we expect the 
above-mentioned quantum phase transition is also quite 
marked. We cannot extrapolate the LSUBn results di- 
rectly, because the phase transition points shift with the 
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order of the LSUBn approximations (Fig. Q). We thus 
find a large statistical deviation of the extrapolated re- 
sults in the region —1.4 < J' < —1.0. Hence, we use 
the minima for M in that region to extrapolate an esti- 
mation of the order parameter. We find that minimum 
[M w 0.05) to be at a value of J' —1.2. The extrapo- 
lated ED results do not agree very well with this result, 
since these give M to be zero at J' w —0.8. However, 
these results are also very poor in that region because 
of the strong influence of the boundary conditions and 
large statistical errors. As a result of these difficulties 
we are not able to decide whether or not quantum fluc- 
tuations and frustration are able to form a disordered 
quantum spin liquid phase (i.e., with M = 0) between 
the Neel state and the helical state for some finite frus- 
trating J' < regime. However, the CCM results suggest 
that there is either no quantum spin liquid phase or that, 
if it does exist, it does so only in a very small region. 

V. SUMMARY 

Using the CCM we have studied the influence of quan- 
tum spin fluctuations on both the ground-state phase 
diagram and the excited states of a spin-half square- 
lattice Heisenberg antiferromagnet with two kinds of 
nearest- neighbour exchange bonds. The phase diagram is 
found to contain a quantum helical phase, a Neel ordered 
phase, and a finite-gap quantum disordered phase. While 
we have clearly a second-order transition from the Neel 
phase to the finite-gap quantum disordered phase, we also 
found indications of a quantum-induced first-order tran- 
sition from the Neel phase to the helical phase, for which 
classically we have a second-order transition. While our 
CCM results were in general in good agreement with the 
ED data, we found the CCM particularly good at de- 
scribing-, the dimerized phase. By contrast, spin wave 
theoryEil fails in that region due to enhanced longitudi- 
nal spin fluctuations. Accurate high-order CCM results 
for the antiferromagnet on the honeycomb lattice were 
also presented. 
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TABLE I. Number of fundamental ground-state configu- 
rations of the LSUBn approximation for the Hamiltonian of 
Eq. ([[]), using a Neel state ($ = 0) and a helical state ($ / 0) 
for the CCM model state, and the number of fundamental ex- 
cited state configurations using the Neel model state only. 



LSUBn 


ground state: $ = 


$ / 


excited state: $ = 


2 


3 


5 


1 


4 


22 


76 


16 


6 


267 


1638 


331 


8 


4986 


42160 


7863 



TABLE II. Ground-state energy per spin, sublattice mag- 
netization, and excitation energy gap for the Hamiltonian of 
Eq. (Q) with J' = 0. This special case is equivalent to the 
honeycomb lattice. We present (extrapolated) LSUBn results 
and extrapolated ED results. 





LSUB2 


LSUB4 


LSUB6 


LSUB8 


extrapolated 


ED 


E/N 
M 
gap 


-0.525 
0.399 
1.182 


-0.540 
0.354 
0.678 


-0.542 
0.334 
0.476 


-0.543 
0.321 


-0.5447 
0.28 
0.02 


-0.543 
0.23 
0.06 
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FIG. 1. Illustration of the classical helical state for the 
square-lattice Heisenberg antiferromagnet of Eq. (Jl]), with 
two kinds of regularly distributed nearest-neighbour exchange 
bonds, J (dashed lines) and J' (solid lines). The spin orien- 
tations at A and B lattice sites are defined by the angles 
8 n = n$ and n = n$ + n, respectively, where n = 0, 1, 2, 
and $ is the characteristic angle of the helical state. The state 
is shown for $ = 7r/12 and n = 0, 1, . . . , 7. 
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FIG. 2. Ground-state magnetic order parameter (Eq. (^)) 
versus J', for the CCM-LSUBn approximation. The results 
are compared (for the Neel region only) with Af s (oo), using 
exact diagonalization (ED) data 

of the antiferromagnetic structure factor, using the ansatz 



(1/7V 2 )V\ (si-Sj) = M s (oo) 2 +constxA 



-1/2 



Note that both extrapolated results fit poorly in a region 
around J' f» —1, and we therefore plot them here as isolated 
points (omitting the solid lines). 
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FIG. 3. The gap A between the lowest-lying excitation 
energy and the ground-state energy versus J' using the 
CCM-LSUBn approximations, in comparison with the extrap- 
olated result of exact diagonalization (ED) (using the ansatz 
A = Aoo + const x TV" 1 ). 
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FIG. 4. Ground-state energy (Eq. (Q) versus J' for 
the extrapolated CCM-LSUBn approximations, in com- 
parison with results of ifmin-wave theory (linear and sec- 
ond-order renormalized) ,H and with the extrapolated re- 
sult of exact diagonalization (ED) data (using the ansatz 
E/N = E^/N + const x N~ s/2 ). 
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FIG. 5. Ground-state structure fac- 

tor S(k) oc ^\ e^ R 3' -R "' k (sj ■ Sj) (i.e., the summation is 
taken over one sublattice) for the Hamiltonian of Eq. (|l|) with 
32 spins, for the quantum and the classical case. The Neel or- 
der (k = (0,0)) becomes unstable against the helical order in 
the classical model for J' < —0.5, but in the quantum model 
the Neel ordering gives way to helical order only for J' < —1.1 
(i.e., the Neel ordering is stable quantum-mechanically in a 
region where it is classically already unstable). 
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FIG. 6. Ground-state energy of the Hamiltonian of Eq. ([j]) 
using CCM-LSUB4 versus the parameter $ of the helical 
CCM model state for certain values of J' in the range 
-1.5 < J' < -1.0. A local minimum of E($) at $ / ap- 
pears for J' < —1.1, which for values of J' < —1.35 becomes 
a global minimum (i.e., at $ = 3>lsub4(J'))j indicating the 
typical scenario of a first-order phase transition. The arrows 
(v) indicate the value of 4> c i for the different values of J'. 
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FIG. 7. The angle "^LSUBn (which minimizes the energy 
-E($), see Eq. (p3[)) versus J', compared with the correspond- 
ing classical result <3> c i (see Eq. (Q)). We find in the quan- 
tum case (LSUBn) a first-order phase transition (e.g., for 
LSUB6 at J' ~ —1.35 where ^LSUBn jumps discontinuously 
from zero to about 0.65). By contrast, in the classical case a 
second-order transition occurs at J' = —1/3. 
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